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The mass and width of the W boson are determined in e + e~ collisions at LEP using 183 pb _1 of 
data recorded at a centre-of-mass energy \fs = 189 GeV with the OPAL detector. The invariant 

Qh| mass distributions from 970 W + W~ — ► qqqq and 1118 W + W~ — ► qq£f£ candidate events are used to 
measure the mass of the W boson, 

M w = 80.451 ± 0.076(stat.) ± 0.049(syst.) GeV. 

"X' 

A direct measurement of the width of the W boson gives Tvv = 2.09 ± 0.18(stat.) ± 0.09(syst.) GeV. 

The results are combined with previous OPAL results from 78 pb" 1 of data recorded with yfs from 

161 to 183 GeV, to obtain: 

M w = 80.432 ± 0.066(stat.) ± 0.045(syst.) GeV, 
T w = 2.04 ± 0.16(stat.) ± 0.09(syst.) GeV. 

The consistency of the direct measurement of Myv with that inferred from other measurements of 
electroweak parameters provides an important test of the Standard Model of electroweak interactions. 
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1 Introduction 



The LEP e + e~ collider at CERN provides an ideal environment for the study of the properties of 
the gauge bosons of the Standard Model (SM) of electroweak interactions [pj. In the first stage of 
its operation, LEP produced e + e~ collisions at centre-of-mass energies, y/s, within a few GeV of the 
Z° resonance, allowing precise measurements of the properties of the Z° boson including its mass and 
fermionic couplings [g,|3[ . In the context of the SM these measurements place constraints on the mass 
of the Higgs boson and provide an indirect determination of the mass of the W boson, M\y- Since 
1996, the LEP collider has operated above the threshold for W + W~ production (LEP2), allowing 
measurements of the trilinear gauge boson couplings Q and a direct measurement of M\y. When 
combined with the direct measurements of the top quark mass at the Tevatron , measurements of 
Mw enable further constraints to be set on the mass of the Higgs boson via electroweak radiative 
corrections H, M. Comparison between the direct measurements of the mass of the W boson and 
the value determined indirectly from data recorded at -y/i as Mi provides an important test of the 
self-consistency of the Standard Model. The direct measurement of Tw further tests the consistency 
of the Standard Model. 

The combination of direct measurements of M\y from LEP2 at yfs ~ 161 — 183 GeV |J8 - 13 1 and 
from hadron colliders [14] currently give il% = 80.419 ± 0.056 GeV [15|. This direct measurement 
is consistent with the indirect value obtained from lower energy data, primarily measurements at 
^/s ~ M z , which give M$ direct = 80.382 ± 0.026 GeV |]. 

The published OPAL measurements of M\y are based on approximately 78 pb _1 of data. This 
paper describes a measurement of the mass of the W boson using a further 183 pb _1 of data recorded 
by OPAL during 1998 at y / s~ 189 GeV. This result is combined with previous OPAL measurements 
to give a direct measurement of the mass of the W boson with a total uncertainty of 79 MeV. It is 
expected that the ultimate LEP precision on M\y will be approximately 30 MeV when all data are 
included and the results of the four LEP experiments are combined [16]. 



2 Data and Monte Carlo Samples 



A detailed description of the OPAL detector can be found in [17]. The data sample used for this anal- 
ysis corresponds to an accepted integrated luminosity, evaluated using small angle Bhabha scattering 
events observed in the forward calorimeters fl8fl , of 183.1±0.2(stat.) ±0.4(syst.) pb _1 . The luminosity 
weighted mean centre-of-mass energy for the data sample is y/s = 188.635 ± 0.040 GeV IT 



2.1 Monte Carlo samples 

A number of Monte Carlo generators are used to simulate the physics processes relevant to the stud- 
ies presented in this paper. All samples include a full simulation of the OPAL detector ^Cj. The 
main physics processes at LEP2 can be broken down into three main categories: four-fermion (4/) 
production, including e + e~ — > W + W~ — > 4/ and e + e~ — > Z°Z° — > 4/, but excluding contributions 
from multi-peripheral diagrams; two-fermion production; and multi-peripheral two-photon mediated 
processes. For the measurement of the W boson mass and width, only the four-fermion processes and 
the two-fermion background process e + e~ — > Z /-/ — > qq play an important role. 



The Koralw program |21], which uses matrix elements calculated with grc4f [22], is used to 
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simulate the production of most four-fermion final states, e + e~ — ► 4/. The main Monte Carlo samples 
are generated at y/s = 188.634 GeV with M\y = 80.33 GeV. The samples were generated using 
the running width scheme for the Breit-Wigner distribution. The four-fermion samples are divided 
into final states which have contributions from processes involving the W-boson propagator and those 
which do not. In the invariant mass region close to M\y, the four-fermion cross-section is dominated 
by doubly resonant W-pair production diagrams (CC03)[|. Additional Koralw four-fermion Monte 
Carlo samples are produced with different centre-of-mass energies and with different values of M\y. 
The four-fermion background from the process e + e~ — > e + e~qq is simulated using the grc4f p^| 
generator. The most important two-fermion background process, e + e~ — ► Z°/7 — ► qq, is simulated 
using Pythia [23], with Herwig [Q used to assess possible systematic uncertainties. 



3 W+W~ Event Selection 



The event selections are described in [|25|] and references therein. The selections are sensitive to the 
leptonic W + W~ — » Ivftvg, semi-leptonic W + W~ — > qq£i^£, and hadronic W + W~ — > qqqq final states. 
Due to the presence of two unobserved prompt neutrinos in the W + W~ — > iviivi final state there is 
little sensitivity to Jl% and the leptonic final state is therefore not used here. 

Semi-leptonic W + W~ — > qqfe^ decays comprise 44% of the total W + W~ cross-section. The event 
selection employs three multivariate relative likelihood discriminants, one for each of the W + W~ — ► 
qqef e , W + W~ — ► qq//i/^, and W + W~ — > qqrv T final states. The W + W~ — ► qqei/ e and W + W~ — > 
qq^f^ channels are characterised by two well-separated hadronic jets, a high-momentum lepton and 
missing momentum due to the prompt neutrino from the leptonic W decay. The signature for the 
W + W~ — > qqrv T channel is similar, with the exception that the r lepton is identified as an isolated, 
low-multiplicity jet typically consisting of one or three tracks. W + W~ — ► qqii'i events are selected 
with an efficiency of 87% and a purity of 91%. The dominant backgrounds are Z°/7 — v qq and 
four-fermion processes such as e + e~ — ► Wef e and e + e~ — » (Z°/7)*(Z°/7)* — > qq£ + ^ _ . 

Hadronic W + W~ — > qqqq decays comprise 46% of the total W + W~ cross-section and are identified 
by requiring four energetic hadronic jets and little or no missing energy. A preselection removes 
approximately 98% of the dominant background process, Z°/7 — > qq. A multivariate relative likelihood 
discriminant is then employed to select the W + W~ — ► qqqq candidates with an efficiency of 87% and 
a purity of 77%. 

After the selections are applied, 1546 W + W~ — ► qqqq and 1246 W + W~ — > qqiise candidate events 
remain, consistent with Standard Model expectations. Not all events are used in the measurements 



described here. As discussed in Section |44j, additional cuts are applied to remove poorly reconstructed 
events and further reduce backgrounds. 



4 Measurement of the Mass and Width of the W Boson 



The measurement of the mass and width of the W boson proceeds in two stages. Firstly the invariant 
masses of the W decay products are reconstructed on an event-by-event basis. Kinematic fits are 
applied to each selected W + W~ — ► qqii^i and W + W~ — ► qqqq event to improve the mass resolution. 
The reconstructed invariant mass spectra are then used to determine Jl% and Tw- 
in this paper, the doubly-resonant W pair production diagrams, i.e. t-channel v e exchange and s-channel Z°/7 
exchange, are referred to as "CC03", following the notation of lid]. 



Fits to the invariant mass spectra to obtain M\v are performed using three different techniques. 
The central results of this paper are obtained using a Monte Carlo reweighting technique |]16| to 
fit the observed mass spectra to obtain J% and I\v- The W mass is also determined using two 
alternative methods, which are used as cross-checks. In the first, an analytic fit to the measured mass 
spectrum uses an unbinned likelihood fit method to determine Jl%. To describe the signal shape, the 
fit uses a parametrisation based on a Breit-Wigner function M . The second method uses a convolution 



technique [12, T^]. The three different fitting techniques have similar expected statistical sensitivities 



and similar estimated systematic uncertainties. 



4.1 Invariant mass reconstruction 



The three methods for extracting M\y use nearly identical procedures to reconstruct the invariant 
mass of the W candidates. The description here applies to the reweighting method. Small variations 



relevant for the alternative analyses are discussed in Section 4.3. 



In previous OPAL measurements of M\y the tracks and clusters in selected W + W~ — * qqqq 
events were grouped into four jets using the Durham algorithm |26|| . For the results presented here, 
W + W~ — ► qqqq events are grouped into either four or five jet topologies depending on the value of 
ln(?/45), where 2/45 is the value of the Durham jet resolution parameter at which the transition from 
5 to 4 jets occurs. Events with ln(y45) > —5.6 are treated as five jets. This separation allows for 
the possibility of hard gluon radiation from one of the quarks and is found to improve the expected 
statistical sensitivity to Mw by approximately 5% compared with treating all events as four jets. The 
improvement in statistical sensitivity comes from events where, under the four jet hypothesis, a gluon 
jet from one W boson is combined with a quark jet from a different W boson. In the semi-leptonic 
decay channels, W + W~ — ► qqtvi, the lepton candidate is removed and the hadronic part of the event 
is reconstructed as two jets. The division of W + W~ — ► qq£u£ events into two and three jet would 
not improve the mass resolution. After the association of tracks and clusters into jets, corrections 
derived from the Monte Carlo simulation are applied to the measured jet momenta to account for 



double counting from particles which deposit energy in more than one sub-detector [27|. 



The invariant masses of the two W bosons can be determined directly from the reconstructed 
momenta of the observed decay products. However, the mass resolution is limited by the relatively 
large uncertainty on the measured energies of the jets, <je/E ~ 20%, rather than by uncertainties 
on the measured jet directions. For this reason the use of a kinematic fit which imposes the four 
constraints of energy and momentum conservation |28J (4-C fit) significantly improves the invariant 
mass resolution. The improvements mainly arise from the energy in event being constrained to the 
well measured centre-of-mass energy. The output of the 4-C fit consists of two reconstructed masses 
per event, one for each W boson in the final state. The resolution of the kinematic fit can be improved 
by imposing a further constraint that the masses of the two reconstructed W boson candidates are 
equal (5-C fit), yielding a single reconstructed mass per event. For W + W~ — ► qqlvi events the three 
unmeasured variables corresponding to the neutrino momentum means that the effective number of 
constraints in the semi-leptonic is two, giving a 2-C fit. For the results presented in this paper, the 
mass reconstruction is performed using the 5-C fit for W + W~ — > qqqq events and the 2-C fit for 
W + W~ — > qqff^ events. 

A common kinematic fitting algorithm is used for W + W~ — > qqef e , W + W~ — » qq/xf^ and 
W + W~ — ► qqqq events. The fitted mass is obtained using an iterative x 2 - m i n i m i sa ti° n procedure 
where the constraints in the kinematic fit are implemented using Lagrange multipliers. The presence 
of initial state radiation (ISR) is neglected. The corrected jet momenta, their associated errors, and 
the measured jet masses are input to the kinematic fit. The errors on the measured jet momenta 



are parameterised by expressions derived from Monte Carlo studies, which are functions of the visible 
energy and polar angle of the jet. The jet masses are fixed to their measured values. This is found to 
improve the mass resolution from the fit compared to treating jets as massless and also to reduce the 
bias in the fitted mass distribution. 

The mass reconstruction for W + W~ — ► qqr^ T events is different. In the W + W~ — - ► (\c\tv t sys- 
tem, the absence of a measurement of the tau lepton energy means that the W mass information is 
determined entirely by the hadronic system. The absence of a measurement of the tau lepton energy 
reduces the effective number of constraints in the kinematic to one. Previously [10] the reconstructed 



mass in W + W~ — ► qqrz/ T events was obtained from the invariant mass of the jet-jet system, scaled by 
the ratio of the beam energy to the sum of the jet energies. For the results presented here an analytic 
formula which reproduces the results of the 1-C kinematic fit is used. 

For W + W~ — » qqef e and W + W~ — ► qqpiVn events the lepton direction is taken to coincide with 
the direction of the track associated with the electron or muon candidate. The energy is estimated 
from the associated electromagnetic calorimeter cluster for electrons and from the momentum of the 
track for muons. Unassociated electromagnetic clusters close to the lepton track, consistent with being 
from final state radiation (FSR), are included in the energy calculation. For each event a single mass is 
determined from a 2-C kinematic fit. In addition to the fitted mass, m iec , the error on the fitted mass, 
cr rec , and the x 2 fit probability, -Pfits are calculated. Events with m rec > 65 GeV and Pfit > 0.001 
are retained. About half of the W + W~ — ► qqriv events selected as either W + W~ — ► qqef e or 
W + W~ — > qq/if^ fail these. The mass information is recovered by treating events with m TCC < 65 GeV 
or PpiT < 0.001 as W + W _ — > qqr^v events. In addition, about 4% of selected W + W~ — ► qqef e and 
W + W~ — ► qq/xf^ events have the identified lepton beyond the effective tracking acceptance of the 
OPAL detector, within 20° of the beam axis, in which case the lepton energy is either poorly measured 
or not measured at all. These events are included in the W + W~ — » qqri/ r channel. The numbers of 
events used in each of the W + W~ — ► qq£ve channels are given in Table |l|. 

The situation for W + W~ — ► qqqq events is complicated by the fact that there are three possible 
assignments of four jets to the two W bosons. For events reconstructed as five jets there are ten possible 
assignments of the jets to the W bosons. Incorrect combinations contain little or no information on 
the mass of the W boson. The incorrect jet-pairings result in a combinatorial background. For each 
W + W~ — > qqqq event, three (or ten) kinematic fits are performed, corresponding to the possible jet- 
pairings. To eliminate poorly reconstructed events and reduce backgrounds, only combinations which 
give a successful 5-C kinematic fit with a resulting x 2 fit probability Pfit > 0.01 and m rcc > 65 GeV 



are considered. A multivariate relative likelihood discriminant similar to that described in [1C] is 
employed to pick out a single combination for each event and reduce combinatorial background. 
Different variables are used for the jet-pairing likelihoods for the four and five jet cases. For four 
jet events, the combination used in the fit is selected on the basis of two variables: the difference 
between the two fitted masses from the 4-C fit and the sum of the di-jet opening angles. In the five 
jet sample four variables are used: the 5-C fit mass, the difference between the two fitted masses from 
the 4-C fit, the minimum opening angle between the jets in the system assigned as W^ — > qqg and 
the cosine of the polar angle of the reconstructed W — ► qqg system. For events reconstructed as four 
(five) jets, the combination corresponding to the largest jet-pairing likelihood is retained provided it 
has a likelihood output exceeding 0.40(0.42). Monte Carlo studies indicate that in 89% (70%) of the 
surviving signal events, the selected combination corresponds to the correct jet-pairing. The number of 
surviving events in the W + W~ — » qqqq channel is given in Table |l|. Figure [y shows the reconstructed 
mass distributions in both four and five jet channels before and after the jet-pairing likelihood cuts. 



For events without ISR the average m rec resolution (as defined in |10||) for the correct jet pairing 
in W + W~ — ► qqqq events is 1.7 GeV. For W + W~ — ► qq^ events the average m rec resolution is 
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2.4 GeV, 2.8 GeV and 3.4 GeV in the W+W" -► qqei/ e , W+W~ -♦ qqfjtu^ and W+W~ -> qqru T 
channels respectively. 



4.2 Extraction of the W mass and width 

The Monte Carlo reweighting technique is used to provide the central results of this paper. The W 
boson mass and width are measured by directly comparing the reconstructed mass spectra in the data 
to Monte Carlo mass spectra corresponding to different values of M\y an d I\y- A likelihood fit is 
then used to extract Jl% and Tw by determining which Monte Carlo spectrum best describes the 
data. The Monte Carlo spectra for arbitrary (My/, I\v) are obtained using the reweighting technique 
described in ||]. In previous OPAL publications, W + W~ Monte Carlo samples generated using only 
CC03 diagrams were reweighted using the ratio of Breit-Wigner functions. For this analysis, Koralw 
e + e~ — ► 4/ reference samples generated with (Jl% = 80.33 GeV, T\y = 2.093 GeV) are reweighted to 
(M(y,r^y) using the ratio of Breit-Wigner functions. The reference samples used for the reweighting 
include only final states which have contributions from diagrams involving the W propagator. This 
procedure is found to give a good approximation to the more exact treatment of using the full four- 
fermion matrix elements, introducing a bias of less than 5 MeV in the W + W~ — ► qqef e channel 
and less than 2 MeV in the other channels. Using the ratio of Breit-Wigner functions rather than 
four-fermion matrix elements results in a much faster fit. 

The mass spectra for background events are taken from Monte Carlo and are assumed to be 
independent of M\y and r\y. The main sources of background are e + e~ — ► Z°/7 — > qq and four- 
fermion processes. The background mass distributions are normalised to the expected number of 
background events. The reweighted signal spectra are then normalised such that the total number of 
signal plus background events corresponds to the observed number of events. This is done separately 
for the W + W~ — ► qqqq, W + W~ — ► qqef e , W + W~ — ► qq/UZ/^ and W + W~ — ► qqrtv channels, with 
the W + W~ — » qqqq channel split into four and five jet topologies. In addition, the W + W~ — > qq^ 
channels are divided into four subsamples according to the error on the reconstructed invariant mass, 
cr rec - This division gives a larger weight to events with reconstructed masses which are known with 
better precision (i.e. small a Tec ) and reduces the expected statistical uncertainty on the fitted W mass 
from the W + W~ — ► qqlvi channels by approximately 5%. In the W + W~ — * qqqq channels, the width 
of the reconstructed mass distribution is dominated by the intrinsic width of the W. Consequently a 
similar subdivision into bins of a rec does not improve the M\y sensitivity. However, jet-pairings which 
give a large jet-pairing likelihood are more likely to be correct and have a better mass resolution 
as shown in Figure |2|. For this reason, the W + W~ — > qqqq events reconstructed as four jets are 
subdivided into four bins of jet-pairing likelihood, resulting in a 7% improvement in the statistical 
sensitivity to J%. No subdivision is performed for the five jet events. 

A binned log-likelihood fit to the m rcc distributions of the data is performed in the range m rcc > 
65 GeV. The log-likelihood function is identical to that used previously |K|. A log-likelihood curve is 
determined separately for each channel. For the W + W~ — ► qqlvi channels, the results are obtained 
by adding the log-likelihood curves separately determined from each channel in each bin of o" rec - For 
the W + W~ — > qqqq channel, the results are obtained by adding the likelihood curves obtained from 
the five jet events to the likelihood curves obtained from the four bins of jet-pairing likelihood for 
events reconstructed as four jets. 



Two types of fit are performed. In the one parameter fit, Tw is constrained by the SM relation [jlq| , 

T w = 3G F M^(l + 2a s /3vr)/(2V27r), (1) 



and only A% is determined. The results of this fit for each channel are given in Table |2| and displayed 
in Figure [|. The combined result is discussed in Section y. In the two parameter fit, both Jl% and 
T\v are determined simultaneously 

In the reweighting method the fitted parameters are expected to be unbiased since any offsets in 
the reconstructed mass introduced in the analysis are implicitly accounted for in the Monte Carlo 
reconstructed mass spectra used in the reweighting procedure. This is verified using several Monte 
Carlo samples generated at various Mw and T\y. In addition, tests using a large ensemble of Monte 
Carlo subsamples, each corresponding to 183 pb~ x and including background contributions, are used 
to verify for each channel separately and for all channels combined, that the measured fit errors 
accurately reflect the root-mean-squared spread of the residual distribution for both the Myy and Tyy 
fits. Since M\v in the Monte Carlo corresponds to the running width definition, so does the fitted mass. 
The expected statistical error on the W mass from the combination of the <\<\ivn and qqqq channels 
is 76 ± 1 MeV, where the weights given to the two channels are determined by the both statistical 
and systematic error contributions. The quoted uncertainty on the expected statistical error is from 
Monte Carlo statistics. 



4.3 Alternative fit methods 
4.3.1 Breit-Wigner fit 



The Breit-Wigner method is analogous to that described in |1C]. It employs an unbinned maximum- 
likelihood fit to the reconstructed mass spectrum using an analytic Breit-Wigner function to describe 
the mass spectrum from e + e~ — > W + W~. Due to initial-state radiation, the reconstructed mass 
spectrum is asymmetric. For the W + W~ — > <\<\lvg channels a relativistic Breit-Wigner function, with 
different widths above and below the peak, gives a satisfactory description of the m rec lineshape. The 
fitting function used is 

S(m rec ) = A m ^ r +(-) (2) 

( m rec ~~ rn o) + m rec +(-) 

where r + (_) is the width assumed for all m rec above (below) the peak centred at itiq. This empirical 
choice of fitting function provides an adequate description for samples up to ten times the integrated 
luminosity of the data. The widths, r + (_), are fixed to values determined from fits to W + W _ signal 
Monte Carlo samples. Different widths are obtained for each W + W~ — > qqiue decay channel. The 
shapes of the background distributions and the background fractions are determined from Monte 
Carlo. The background fractions are held constant in the fit. The fit is restricted to the range 
70 < m rec < 88 GeV. 

In the W + W~ — ► qqqq channel events are divided into four and five jet samples which are fit- 
ted separately. The division is made on the basis of the 5-C kinematic fit probability. Events are 
reconstructed as five jets if any of the ten possible jet-pairings in the five jet assignment gives a fit 
probability which is greater than twice that of highest probability of the three jet-pairings for the 
four jet hypothesis. For both four jet and five jet samples, the fitting function S(m Tec ) is multiplied 
by G{m rec ) = exp[— (mo — m rec ) 2 /2a 2 ]. This empirical choice provides a good description of the re- 
constructed mass spectra in the fit range for W + W~ — ► qqqq events. The value of a is determined 
from Monte Carlo. For the four jet W + W~ — ► qqqq sample either one or two jet-pairings are used 
following the procedure described in ||. In the five jet sample, a jet-pairing likelihood is used. Here 



the variables used in the likelihood are those used for the reweighting fit described in Section |4.l| , with 
the exception that the 5-C fit mass is not included. 



In contrast to the procedure employed for the reweighting method, the W + W~ — ► qqtvi events 
are not divided into subsamples according to cr rec nor are the W + W~ — > qqqq events divided into 
subsamples according to the jet-pairing likelihood. However, W + W~ — > qqTv T events are sub-divided 
into fully-leptonic and semi-leptonic decays of the tau lepton. 

The fitted mass, mo, must be corrected for offsets not accounted for in the fit, e.g. from initial-state 
radiation and event selection. The correction is determined using fully simulated Monte Carlo samples 
generated at different values of A% with the expected background contributions included and is found 
to depend linearly on M\y- The results from the W + W~ — > qqqq and W + W~ — » qqlh>i channels, after 
correction, are given in Table y. The expected statistical error on the combined W mass measurement 
is 78 ± 2 MeV. 



4.3.2 Convolution fit 

The convolution method [12,16] attempts to exploit all available information by constructing a likeli- 
hood curve for each selected event. The likelihood is calculated using a functional form 

C(M w ,m rec ) =p s V s (M w ,m Tec ), 

where p s is the probability of a candidate event being a true signal event and V s is the probability 
function for M\y given the observed reconstructed mass, m rec . The function V s (M\M,m Tec ) is defined 
as 

V s {M w ,m rec ) = BW(M w ,m,s') (g>ISR(s,s') ®R(m,m rec ), 

where BW(A%,m, s') is the relativistic Breit-Wigner function for producing off-shell W bosons with 
mass m for a W mass of M\v including the effects of phase space. T\y is fixed relative to M\v using 
Equation [l]. The radiator function, ISR(s, s'), is used to account for the effects of ISR (on a statistical 
basis) . The inclusion of this term accounts for the fact that in the kinematic fit the observed energy is 
constrained to y/s, where in reality the effective centre-of-mass energy, yV, depends on the amount of 
ISR. Finally, R(m, m rcc ) is the probability density function relating the experimentally reconstructed 
mass, m rec , to the average of the two true masses, m, of the off-shell W bosons in the event. 

In the W + W~ — > qqivi channels, the full error information on the fitted mass of each event is used 
in an unbinned maximum likelihood fit. In a previous OPAL publication |p!o[ , the event probability 
function, R(m, m rcc ), was assumed to be Gaussian with its central value and standard deviation 
taken from a 2-C kinematic fit. In this analysis, the 2-C fit is modified, replacing the equal W mass 
constraint by one in which the masses are individually constrained to a given W boson mass, m. The 
X 2 distribution from this fit is converted into the event probability density function, R(m, m TCC ), which 
is non-Gaussian. 

In the W + W~ — > qqqq channel, all events are forced into a five jet topology. A relative likelihood, 
constructed from the difference between the reconstructed masses in a 4-C kinematic fit, is used to 
select on average about three jet-pairings per event thereby reducing significantly the combinatorial 
background. In the Monte Carlo 92% of selected events the correct jet-pairing combination is used in 
the fit. The convolution is performed in two dimensions, namely the reconstructed masses of the two 
W bosons, with joint probability density functions corresponding to each jet-pairing obtained from 
the 4-C kinematic fit. The probability density functions for the different combinations are added to 
form an event probability density function [p 



In both channels the log-likelihood curves from each selected event are summed to yield a single 
curve from which a fitted mass is determined. The fitted mass is corrected for offsets not accounted for 



in the fit in the manner described for the Breit-Wigner fit. The results from the W + W~ — ► qqqq and 
W + W~ — ► qq£f£ channels, after all corrections, are given in Table y. The expected statistical error on 
the combined W mass measurement is 72 ±2 MeV. The total combined expected error, statistical and 
systematic, from the convolution method is slightly smaller than that obtained from the reweighting fit. 
However, at present the reweighting fit is the only method in which a simultaneous fit for (M\v, T\v) is 
implemented. For this reason, and for consistency with previous OPAL publications, the reweighting 
method is retained for the results given in this publication. 



5 Systematic Uncertainties 

The systematic uncertainties for the measurement of Mvv using the reweighting method are estimated 
as described below and summarised in Table ||. The contributions from each of the sources are added 
in quadrature to yield the total systematic uncertainty. All contributions were evaluated for both the 
Mvv an d Tw determinations. For the alternative analyses, the systematics are estimated similarly and 
yield comparable results. The determination of the individual systematic uncertainties is described 
below. 



5.1 Beam energy 

The average LEP beam energy is currently known with a precision of 20 MeV [|l9| which leads to a 
systematic uncertainty of 16 MeV on il% and 2 MeV on Tw- The RMS spread in the LEP centre-of- 



mass energy of 237 ± 12 MeV [19| and the average difference between the electron and positron beam 
energies results in negligible impact on the measurement of Myv- The Monte Carlo does not include 
the intrinsic beam energy spread. This introduces a bias of +10 MeV to the value of T\v obtained 
from the reweighting fit. A —10 MeV correction is applied and the size of the correction taken as an 
additional systematic uncertainty on the measurement I\v. 



5.2 Initial state radiation 

The main Koralw Monte Carlo samples include an 0(a 3 ) treatment of initial state radiation (ISR). 
The systematic error associated with ISR is estimated by reweighting the generated Koralw events 
to correspond to an 0(a 2 ) or 0(a) treatment of ISR using the matrix elements calculated in Koralw. 
The resulting shifts from the comparison of 0(a 3 ) to 0(a) are less than 2 MeV and are neglected. 



5.3 Full 0(ol) corrections 

In the process e + e~ — > 4/7 ISR diagrams dominate. The accuracy of the treatment of ISR in the 
Koralw generator is sufficient for this analysis. However the Koralw generator does not include 
all 0{a) effects. For example, interference between ISR and FSR graphs and photon radiation from 
the W bosons are not implemented. Recent calculations using the double-pole approximation indicate 
that possible mass biases of order 10 MeV could arise |p9[j . However, these estimates were based 
on generator level studies which compared mass distributions from the full O{o) treatment with the 
corresponding Born level treatment which does not include any final state radiation. To obtain a 
proper estimate of the potential mass bias when applied to the experimental measurement it would 
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be necessary to use the full OPAL detector simulation and mass reconstruction procedure. Currently 
this is not possible since the necessary Monte Carlo programs are not yet available. For this paper 
possible systematic biases from the full 0(a) treatment of photon radiation are not included since the 
size of the effect is significantly smaller than the current statistical precision of the M\y measurement. 



5.4 Hadronisation in W 1 " 1 — ► qq decays 

The Koralw Monte Carlo samples taken as references in the extraction of Mw use the Lund string 
model as implemented in Jetset for the simulation of the hadronisation^ of W — ► qq decays. The 
Jetset model was tuned to OPAL hadronic data recorded at the Z° resonance [pQl . The model was 
tuned to describe event shape variables and inclusive particle production rates. To assess a possible 
systematic error from the uncertainties in the tuning of the Jetset model, the Jetset parameters 



aq, b, Aqcd and Qq are each varied in turn by ±1<t about their tuned values at the Z° [30]. The 
maximum resulting biases to the fitted values of M\y from these high statistics samples are 30 MeV 
and 20 MeV in the qqivt and qqqq channels respectively. In addition, an earlier OPAL tune of the 
Jetset model is used. This version was tuned to describe global event shape variables alone |3l| . The 
resulting differences in the fitted value of _M\y using the earlier OPAL tune are (—33 ± 14) MeV and 
(—30 ± 12) MeV in the W + W~ — > qqqq and W + W~ — > qqivt channels respectively. 

As a further test of hadronisation uncertainties, the Jetset string model is replaced by the Herwig 
cluster model. For this comparison, two samples were generated using the same W + W~ final states and 
differing only in the hadronisation modelling. The resulting fitted mass values from the reweighting 
method from Jetset and Herwig are compared. Differences (Jetset-Herwig) of (—18 ± 14) MeV 
and (—11 ± 12) MeV are obtained for the W + W~ — ► qqlvi and W + W~ — ► qqqq channels. A similar 
comparison is made using the Ariadne colour dipole model (3^]. Differences (Jetset-Ariadne) of 
(+15 ± 14) MeV and (+6 ± 12) MeV are obtained for the W + W~ -+ qqtv t and W+W~ -► qqqq 
channels. 

On the basis of these comparisons systematic uncertainties of 30 MeV are assigned to the W mass 
measurements from both W + W~ — ► qqlvg and W + W~ — ► qqqq. The hadronisation uncertainty is 
taken to be fully correlated between the two channels. For the measurement of T\v the corresponding 
systematic errors in the W + W~ — > qqlvi and W + W~ — » qqqq channels are 25 MeV and 60 MeV 
respectively. 



5.5 Detector calibration and simulation 

The effects of detector mis-calibrations and deficiencies in the Monte Carlo are investigated by varying 
the jet and lepton energy scales over reasonable ranges. The ranges used for the systematic variations 
are dependent on polar angle and are determined from detailed comparisons of data and Monte Carlo 
utilising both data recorded at 189 GeV and 3.1 pb~ of data collected at y/s ss Mz during 1998. 
From e + e~ — ► e + e~ events the electro-magnetic calorimeter energy scale is known to 0.3%. Similarly 
from e + e~ — ► Z° — ► /i + /-i _ events, the uncertainty on momentum scale from the tracking detectors 
is determined to be 0.3%. These scale uncertainties dominate the detector related systematics in 
the W + W~ — ► qqef e and W + W~ — > qqfiv^ channels, and are 14 MeV and 16 MeV respectively. 
Uncertainties on energy (momentum) resolution and angular resolution are also evaluated but have 
much smaller impact. Uncertainties on the jet energy scale are determined from e + e~ — ► Z° — ► qq 

2 In this paper hadronisation refers to the process whereby the quarks produced in W decays produce the hadrons 
observed in the detector including the development of the parton shower. 
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events to be 0.2%— 1.0% depending on the polar angle. The dominant effect of the jet energy scale 
uncertainty is a 6 MeV uncertainty in the W + W~ — ► qqr^ T channel. In the W + W~ — ► qqqq channel 
the sensitivity to the overall jet energy scale is greatly reduced due to the kinematic fit. Here the 
dominant effect is from the global uncertainty in the jet polar angle scale, i.e. the relative length-to- 
width scale of the detector, resulting in a 3 MeV uncertainty on the W mass. The linearity of the 
jet energy scale is checked using e + e~ — ► Z° — ► qq events with a clear three-jet topology. For the 
measurement of Tw the main systematic error arises from uncertainties in the jet energy scale and 
resolution and is 47 MeV in both W + W~ — > qq^ and W + W~ — > qqqq channels. 



5.6 Four-fermion modelling 

Possible systematic effects associated with the modelling of the four-fermion final state, including 
interference between W + W _ diagrams and other four-fermion processes, are investigated by compar- 
ing the grc4f and the Excalibur [33] generators. The comparison is performed using the ratio of 



the Excalibur to grc4f four-fermion matrix elements to reweight the Koralw Monte Carlo sample. 
This procedure avoids introducing an additional statistical error associated with generating additional 
Monte Carlo. The resulting shift in the reconstructed mass is 1 ± 1 MeV in the W + W~ — > qq£v( chan- 
nel and 0±1 MeV in the W + W~ — > qqqq channel. Similar results are obtained when the Excalibur 
Monte Carlo is reweighted using the ratio of the grc4f to Excalibur four-fermion matrix elements. 
Consequently no additional systematic uncertainty is assigned. The corresponding systematic uncer- 
tainties on the measurement of Tyy, of order 30 MeV, are determined by the statistical precision of the 
comparisons. In the W + W~ — ► qqef e channel an additional 5 MeV uncertainty is assigned, associated 
with using the Breit-Wigner functions in the reweighting fit. 



5.7 Background treatment 

Uncertainties in the modelling of background processes mainly arise from the simulation of the 
Z°/7 — ► qq background in the W + W~ — ► qqri/ r and W + W _ — > qqqq channels. A number of sys- 
tematic checks of the normalisation and shape of the Z°/j — > qq background are performed. The 
background normalisation in the reweighting fit is varied by one standard deviation, using the uncer- 



tainties evaluated in [25]. The Herwig cluster model is used as an alternative to the Lund string 
model used in Pythia. In addition, for the W + W~ — > qqqq channel, data taken at y/s ~ M%, scaled 
by (189 GeV/Mz), are also substituted for the Z°/7 — ► qq background. The sum in quadrature of 
the above systematic shifts leads to the assignment of a background systematic error of 6 MeV in 
the W + W~ — > qqqq channel and 8 MeV in the W + W~ — > qqtvi channel. The systematic error on 
Tw due to uncertainties in background modelling is 25 MeV and 41 MeV in the W + W~ — » q_q£vi 
and W + W~ — » qqqq channels respectively arising from differences between Jetset and Herwig and 
uncertainties on the accepted background cross section. 



5.8 Monte Carlo statistics 

The finite statistics of the Monte Carlo samples used in the reweighting procedure result in an ad- 
ditional uncertainty of 10 MeV on the W mass determined separately in the W + W~ — > qqqq and 
W + W~ — ► qqlve channels. The corresponding systematic uncertainties for the W width measurement 
are 51 MeV. 
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5.9 Final state interactions (FSI) 

A significant bias to the apparent W mass measured in the W + W~ — ► qqqq channel could arise if the 
hadronisation of the two W bosons is not independent [16|, 34|. At LEP2 energies, the decay length for 
W bosons is about 0.1 fm, significantly less than the typical hadronisation scale, approximately 1 fm. 
The hadronisation of the two W bosons in the W + W~ — ► qqqq channel may have space-time overlap 
allowing possible final state interactions between the decay products of different W bosons. Two 
possible sources of final state interactions, both leading to non-independent hadronisation, have been 
widely considered: Colour Reconnection (CR) and Bose-Einstein Correlations (BEC). Whether these 
effects play a significant role at LEP2 has yet to be experimentally determined. As a result, to assign a 
systematic error related to possible final state interactions it is necessary to rely on phenomenological 
models. The strategy adopted here is to consider only those models for CR and BEC which are 
consistent with both LEP1 and LEP2 data. 

The CR systematic error refers to possible biases in the W mass measurement arising from QCD 
reconnection effects in the non-perturbative phase. Reconnection effects in the perturbative phase, 
i.e. multiple gluon exchange between quarks from different W bosons, have been shown to be small 
[Pq|. To investigate the systematic biases originating from non-perturbative CR effects, the models 
implemented in the Pythia and Ariadne [32] Monte Carlo generators are studied. In both cases the 



results from the mass fit to a Monte Carlo sample including CR effects are compared to the results 
obtained from the corresponding sample without colour reconnection. The results are summarised 
in Table || The largest biases are seen in the Ariadne 2 and Ariadne 3 models]]. However, these 
models have been strongly disfavoured by studies of three jet events in LEP1 Z° data [37]. For this 



reason the current Ariadne implementation of CR is not used to assign a systematic uncertainty. 
Of the remaining models considered, the SK I model [35] produces the largest bias. The SK Monte 
Carlo samples use the Jetset tune of pH], where the parton shower cut-off, which directly influences 



the behaviour of the CR model [35], was set to 1.9 GeV. The strength parameter ki (also called p) 
of the SK I model was set to 0.9, corresponding to 35% of events being reconnected. A systematic 
uncertainty of 66 MeV on M\y is assigned. The corresponding uncertainty on T\y is 68 MeV. 

As is the case in Z° decays, BEC between like-sign charged pions are observed in W^ decays at 
LEP2 p8|. These correlations are not implemented in the standard Monte Carlo programs used in 
this analysis. Correlations between pions from the same W do not affect significantly the global jet 
properties and are unlikely to bias the Jl% measurement. If BEC exist between pions from different 
W bosons in W + W~ — > qqqq events, significant biases to the reconstructed value of il% may arise. To 
investigate possible systematic biases from BEC, a sample of W + W~ — ► qqqq events is produced using 
the Pythia Monte Carlo generator which includes BEC as described in |3JJ. For the evaluation of the 
systematic error the BEC model BE32 of [39] with A = 1.0 and R = 0.42 fm is used. The difference in 



the fitted W mass between a sample generated with BEC between different W bosons and a sample 
generated with BEC only between particles from the same W is +67 ± 14 MeV (different-same). The 
resulting systematic uncertainties on il% and Tyj from BEC in the W + W~ — ► qqqq channel are 
67 MeV and 39 MeV respectively. 



5.10 Fit procedure 

A comparison of the results from the reweighting method with those from the Breit-Wigner and con- 
volution fits, summarised in Table ||, is used to test for residual biases in the reweighting method. The 



Following the notation used in [p6|, |37[ . 
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expected RMS of the differences of the results from the three methods are determined using an ensem- 
ble of Monte Carlo subsamples including background, each corresponding to an integrated luminosity 
of 183 pb _1 . For each subsample the difference in the fitted Jl% determined using the reweighting 
method and that determined using each of the other methods is calculated for the W + W~ — ► qqqq 
and W + W _ — ► qqivi channels. These differences have mean values consistent with zero and RMS 
values of approximately 54 MeV and 41 MeV (62 MeV and 84 MeV) when comparing the reweighting 
fits to the Breit-Wigner (convolution) fits in the W + W~ — ► qqi^e and W + W~ — > qqqq samples, re- 
spectively. Since the alternative analyses yield results consistent with those obtained using the default 
reweighting analysis, no additional systematic error is assigned. 



6 Results 



For the reweighting method described in Section 4.2 , the results of a simultaneous fit to M\y and Tvv 
from the combined W + W~ — > qqqq and W + W~ — > qqtvi, event samples are 

M w = 80.451 ± 0.076(stat.) ± 0.050(syst.) GeV, 
T w = 2.09 ±0.18(stat.)±0.09(syst.) GeV. 

The correlation coefficient between M\y and Tw is 0.04. For this fit, the central values are determined 
by adding the log-likelihood curves from the W + W~ — > qqqq and W + W _ — » qq£v£ channels including 
the effects of systematic uncertainties. The likelihood contours for this fit are displayed in Figure ||. 

A one parameter fit for the mass is performed by constraining the width using Equation || to 
give M w = 80.402 ± 0.104(stat.) ± 0.101(syst.) GeV in the W+W" -> qqqq channel, and M w = 
80.478 ± 0.104(stat. ) ±0.038(syst.) GeV in the W+W~ -> qqiut channel. The combined result, taking 
into account the correlated systematics between the W + W~ — > qqqq and W + W~ — » qqlvi channels, 
is 

M w = 80.451 ± 0.076(stat.) ± 0.049(syst.) GeV. 

Due to the larger systematic errors in the W + W~ — ► qqqq channel, mainly from FSI, it carries a 
reduced weight of 0.36 in the combination. The combined W + W~ — ► qqqq and W + W~ — ► qqlu^ 
results from the alternative analyses, after all corrections, are for the Breit-Wigner fit, Mw = 
80.436 ± 0.078(stat.) ± 0.050(syst.) GeV, and for the convolution fit, M w = 80.378 ± 0.073(stat.) ± 
0.051(syst.) GeV. 



6.1 Combination with previous data 

The results presented in this paper, based on 183 pb _1 of data, are combined with previous OPAL 
measurements of Mvv from direct reconstruction at y^ ss 172 GeV |9| and at y^ « 183 GeV [10]. The 
results in the W + W~ — * qqlvi and W + W~ — > qqqq channels are combined independently taking into 
account systematic uncertainties which are correlated between data from the different centre-of-mass 
energies, giving: 



M W (W + W~ -> qq£v e ) = 80.441 ± 0.086(stat.) ± 0.034(syst.) ± 0.017(lep) GeV, 
M W (W+W" -» qqqq) = 80.409 ± 0.093(stat.) ± 0.034(syst.) ± 0.016(lep) ± 0.094(fsi) GeV. 

where the uncertainties from the LEP beam energy (lep) and from final state interactions (fsi) are 
quoted separately. In the combination the systematic uncertainties for the colour reconnection and 
Bose-Einstein correlations of the previous data were set to those presented in this paper. 
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The difference between the fitted Jl% in the W + W — > qqqq and W + W — > qqlvi channels is 
AM W = (M{f m - Af$^) = -0.032 ± 0.127(stat.) ± 0.047(syst.) GeV. A significant non-zero value 
for AM\y could indicate that FSI effects are biasing the value of M^ determined from W + W~ — » 
qqqq events. The systematic error on the quoted value of AMw does not include contributions 
from CR/BEC effects. For the AM\y measurement the hadronization systematic uncertainties in the 
W + W~ — » qqqq and W + W~ — > qqlvt channels are taken to be uncorrelated; a conservative choice in 
the absence of a good estimate of the true correlation coefficient. 

The Jl% results from direct reconstruction in the W + W~ — > <\<\lvg and W + W~ — > qqqq chan- 
nels are combined with the OPAL measurement of Jl% from the W + W~ production cross-section at 
y/s ~ 161 GeV pi. The combination is made assuming that the mass measurements from direct recon- 
struction and from the threshold cross-section are uncorrelated, apart from the uncertainty associated 
with the LEP beam energy, which is taken to be fully correlated. The direct reconstruction measure- 
ments from the W + W~ — > qqlvi and W + W~ — > qqqq channels from three different centre-of-mass 
energies are combined using the full covariance matrix. The combined result is 

M w = 80.432 ± 0.066(stat.) ± 0.032(syst.) ± 0.028(fsi) ± 0.017(lep) GeV. 
The measurement of T\y presented in this paper is combined with the previous OPAL measure- 



ments from direct reconstruction [H, 10 1 to give 



T w = 2.04 ± 0.16(stat.) ± 0.09(syst.) GeV. 

For the earlier data, with relatively low statistics, the measured statistical error is correlated with the 
measured width. Consequently, to avoid biasing the combination, the separate width measurements 
are weighted using the expected statistical error rather than the measured statistical error. 



7 Summary 

Using the 183 pb _1 of data recorded by the OPAL detector at a mean centre-of-mass energy of 
approximately 189 GeV, a total of 2088 W + W~ — ► qqqq and W + W~ — ► qq^ candidate events are 
used in a fit constraining T\y using Equation [l] to obtain a direct measurement of the W boson mass, 
Mw = 80.451 ± 0.076(stat.) ± 0.049(syst.) GeV, while a second fit is used to determine directly the 
width of the W boson, T w = 2.09 ± 0.18(stat.) ± 0.09(syst.) GeV. 

The results described in this paper are combined with the previous OPAL results from data 



recorded at y/s ~ 161 — 183 GeV flg-jLOj. From this combined data sample the W boson mass is 
determined to be 

M w = 80.432 ± 0.066(stat.) ± 0.045(syst.) GeV. 

This result is consistent with both the other direct measurements and the indirect value inferred 
from fits to electroweak data ||. It should be noted that the systematic error is dominated by the 
uncertainties arising from hadronisation and final state interactions. For future LEP combinations 
these uncertainties will need to be reduced to benefit from the ultimate LEP combined statistical 
uncertainty of about 25 MeV. 

The result for the W boson width is combined with the previous OPAL result from data recorded 
at y/s « 172 - 183 GeV to obtain 

T w = 2.04 ± 0.16(stat.) ± 0.09(syst.) GeV. 
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This measurement of T\v is in good agreement with the Standard Model expectation and with other 



direct measurements at LEP [12, 13, 5J and at the Tevatron [41|. 
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Channel 


Observed 


Expected 


Purity 


W+W" - 


> qqqq (4-jet) 


701 


694 


83% 


w+w- - 


> qqqq (5-jet) 


269 


277 


91% 


w+w-^ 


qqef e 


350 


376 


98% 


w+w-^ 


qq/^ 


365 


373 


99% 


W+W"^ 


qqrv T 


403 


403 


89% 


Combined 


2088 


2123 


91% 



Table 1: Numbers of events used in the W mass and width determination for each channel 



and all channels combined. Only events surviving the cuts described in Section 4.1 are included. 
The W+W~ — > qqr^ r numbers include some events initially selected as either W+W - — > qqef or 
W+W~ — > qq/ity 



Channel 


Measured M w /GeV 


Expected error/GeV 


W+W" - 


->• qqei/ e 


80.375 ±0.175 


0.164 


w + w- - 


■> qq/u^ 


80.513 ±0.163 


0.168 


w+w- - 


-> qqr^ r 


80.594 ±0.227 


0.220 


w+w- - 


-> qqqq(4-jet) 


80.424 ±0.114 


0.112 


w+w- - 


■* qqqq(5-jet) 


80.290 ±0.257 


0.230 


w+w- - 


->• qq£u£ 


80.478 ± 0.104 


0.104 


W+W" - 


+ qqqq 


80.402 ±0.104 


0.100 



Table 2: Results using the reweighting method for the fit from 183 pb _1 of data taken at y^ ~ 
189 GeV for each of the channels separately and for the combined W+W - — > qqli^i channel. Equation 
is used to constrain T\v to M\y- The expected errors are estimated using an ensemble of Monte 
Carlo subsamples. The errors obtained in the data are consistent with the RMS spread of the errors 
obtained from the Monte Carlo subsamples. 



Channel 



W+W" 
W+W" 



qqqq 



Reweighting 
Measured M w /GeV 



80.48 ±0.10 ±0.04 
80.40 ±0.10 ±0.10 



Breit-Wigner fit 
Measured M w /GeV 



80.47 ±0.11 ±0.04 
80.38 ±0.11 ±0.10 



Convolution fit 
Measured M w /GeV 



80.38 ±0.11 ±0.04 
80.37 ± 0.09 ± 0.09 



Table 3: Comparison of M\y fit results obtained using the three different fitting techniques for the 
W+W~ — ► qq^ and W+W~ — ► qqqq channels separately. In each case the first uncertainty is 
statistical and the second systematic. The observed differences are compatible with the expected 
RMS differences between the results from the different fits. 
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Systematic errors 


M w 


r 


w 


(MeV) 


qqqq 


qq^£ 


qqqq 


qq^£ 


Beam Energy 


16 


16 


10 


10 


Hadronisation 


30 


30 


60 


25 


Jet energy/Resolution 


4 


5 


47 


47 


Lepton energy/Resolution 


- 


9 


- 


9 


Four Fermion 


— 


— 


33 


29 


Background 


6 


8 


41 


25 


MC statistics 


10 


10 


51 


51 


Sub-total 


36 


38 


106 


84 


Colour Reconnection 


66 


— 


68 


— 


Bose-Einstein Correlations 


67 


- 


39 


- 


Total systematic error 


101 


38 


132 


84 



Table 4: Summary of the systematic uncertainties for the fit results. For the fits to determine Mw, 
T\v is constrained to its Standard Model relation. The uncertainties are given separately for fits to the 
W + W~ — > qqqq and W + W~ — ► qqii^i samples. Uncertainties from ISR give negligible contributions 
to both M\v and T\y- The uncertainty related to the Monte Carlo models of the four- fermion process 
give is negligible for the measurement of Mw- 



CR Model 


Mass shift (MeV) 


Width shift (MeV) 


SKI 


+66 ±8 


+68 ± 22 


SKII 


+3 ±8 


+20 ± 22 


sKir 


+10 ±8 


+53 + 21 


AR2 


+85 + 8 


+128 + 22 


AR3 


+140 ± 10 


+309 ± 21 



Table 5: The predicted shifts in the fitted mass and width from the W + W~ — ► qqqq channel for 
various colour-reconnection (CR) models. The numbers correspond to the fitted mass (width) using 
the reweighting method in the reconnected sample minus the fitted mass (width) for the corresponding 
sample without CR. The errors are statistical only. 
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Figure 1: The reconstructed mass, m TCC , distributions for W + W~ — ► qqqq events divided into four 
and five jet samples. The top two plots show the distributions before the jet-pairing likelihood (JPLH) 
cut. The lower two plots show the mass distributions which are used to determine Mw, i.e. after 
the cut on JPLH. The points correspond to the OPAL data and the histograms to the Monte Carlo 
predictions. The contribution from the non-WW background is shown as the cross-hatched histogram 
and the combinatorial background is indicated by the singly-hatched histogram. 
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Figure 2: The reconstructed mass distributions for the four jet W + W~ — > qqqq sample divided into 
four bins of the output of the jet-pairing likelihood (JPLH) discriminant. The points correspond to the 
OPAL data and the histograms to the Monte Carlo predictions. The contribution from the non-WW 
background is shown as the cross-hatched histogram and the combinatorial background is indicated 
by the singly-hatched histogram. 
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Figure 3: The reconstructed invariant mass distribution for the W + W~ — ► qqqq, W + W~ — ► qqef e) 
W + W~ — > qq/Ui/^ and W + W _ — v qqrz> r samples. The points correspond to the OPAL data and the 
histograms to the reweighted Monte Carlo spectra corresponding to the fitted masses. The non-WW 
background contribution is indicated by the hatched histogram. 
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Figure 4: The 39% and 86% contour levels of the two parameter fit using the reweighting method 
including systematic contributions. The Standard Model relation between Myy and r\y, given in 
Equation |], is shown by the dashed line. 
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